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We report the magnetic-field induced Bose glass— BEC transition of magnons in Tli-^iK^iCuCls 
and its critical behavior investigated through specific heat and ESR measurements. The field de- 
pendence of the BEC transition temperature Tn can be described by the power law [H — He] oc T^ 
near the quantum critical point ~3.5 T. The critical exponent cj) tends to reach a value smaller 
than 1/2 with decreasing fitting window in contrast with </>—>■ 3/2 for the standard BEC in pure 
system. At sufficiently low temperatures, the ESR line shape for H ~ He is intermediate between 
Gaussian and Lorentzian. This indicates the localization of magnons for H <Hc at T = 0. 

PACS numbers: 72.15.Rn, 75.10.Jm, 75.40.Cx, 76.30.-v 



Heisenberg antiferromagnets composed of spin dimers 
often have gapped ground states at zero magnetic field 
and undergo quantuniphase transitions to ordered states 
in magnetic fields [ll, Q • The field- induced magnetic or- 
dering is typical of the quantum phase transition and 
has been actively studied both theoretically and experi- 
mentally i-IIl.' These studies have established that the 
coupled antiferromagnetic spin-dimer system in a mag- 
netic field can be mapped onto a system of interacting 
lattice bosons, where the bosons are spin triplets with 
Sz = l called magnons or triplons, and that the field- 
induced magnetic ordering is best described as the Bose- 
Einstein condensation (BEC) of magnons. However, for 
the case with exchange disorder, the nature of the ground 
state in a magnetic field and the critical behavior of the 
field-induced magnetic ordering are not sufficiently un- 
derstood, although there have been some studies on this 
problem 

In a magnetic field H comparable to the gap, the ef- 
fective Hamiltonian of magnons is expressed as 

H = ^ (J,; - .g/iB-ff) ajoi 

i 

+ ^v^^Uj + \^Y1 Utjala]aja^. (1) 

i j i j 

The first, second and third terms denote the local po- 
tential, hopping and interaction of magnons, respec- 
tively. The intradimer exchange interaction Ji on dimer 
site i corresponds to the local potential of magnons, 
Vi = Ji—giiBH. The ground state and the quantum 
phase transition for lattice bosons in a random potential 
were investigated theoretically by Fisher et al. [2^ . They 
argued that a new phase called Bose glass (BG) emerges 
as a ground state in addition to BEC and Mott insulat- 
ing (MI) phases, which correspond to the ordered phase 
and gapped phase in the magnetic system, respectively. 
Bosons are localized in the BG phase because of random- 
ness, but there is no gap; thus the compressibility is finite. 
Fisher et al. showed that the BEC transition occurs only 
from the BG phase, and that near T = 0, the relation be- 
tween transition temperature Tc and boson density p is 



expressed as Tc ^ [ps(0)]^, ps{0) ^{p- Pc)'' , where pc is 
the critical density at which the BEC transition occurs 
and Ps(0) is the condensate density at T = 0. Exponents 
y and C are given by y = 3/4 and C > 8/3 for three dimen- 
sions. The critical behavior is different from that of the 
standard BEC case without randomness, for which these 
exponents are y — 2/3 and C = 1- 

Recent theory has demonstrated the emergence of the 
BG phase in a disordered quantum magnet (13, [3 • The 
BG phase has also been studied in other disordered quan- 
tum systems, such as ^He adsorbed on porous Vycor glass 
[2^ [23j . amorphous superconductors jU and trapped 
cold atoms ^5j. From the correspondence between the 
boson density and the magnetization in the coupled spin- 
dimer system, a relation (p — p^.) cc {H — H^) is obtained, 
where He is the critical magnetic field of the magnon 
BEC. Hence, for the coupled spin-dimer system with ex- 
change disorder, the transition field H-[<i{T) near T — 
should be expressed by the power law 

[HN(r)-ffc] cxT^ (2) 

with the critical exponent (j)<l/2. Consequently, the 
low-temperature phase boundary should be tangential 
to the field axis at T~0, as shown in Fig. [T{a). This 
phase boundary behavior is qualitatively different from 
that of pure system, for which = 3/2, and thus, the 
phase boundary is perpendicular to the field axis. Hb 
in Fig. [Ija) is the critical field for the MI-BG transition. 
The BG phase exists between Hb and He- 

To investigate the magnon BEC under the influence 
of the localization, we performed specific heat and ESR 
measurements on Tli_j;Ka;CuCl3. The parent com- 
pounds TlCuCls and KCuCla have the same crystal 
structure composed of the chemical dimer CU2CI6, in 
which Cu"^ ions have spin- 1/2. Their magnetic ground 
states are spin singlets with excitation gaps A/fce of 
7.5 K and 31 K, respectively [2^. The gaps originate 
from the strong antiferromagnetic exchange interaction 
between spins in the chemical dimer. The neighboring 
spin dimers couple antiferromagnetically in three dimen- 
sions. The intradimer exchange interaction was evaluated 
to be J/fcB = 65.9 K and 50.4 K for TlCuClg and KCuClg, 
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FIG. 1: (a) Schematic low-temperature phase boundaries expected for Tli-^^Ki^CuCls. Sohd and dashed hnes denote the 
boundaries for x ^ Q and x = 0, respectively. For simplification, the critical fields He for both the cases are plotted at the same 
position, (b) Field scans of the specific heat in Tli-ajKjjCuCla with a; = 0.36 measured at various temperatures. Arrows denote 
the transition field H]si{T). (c) Magnetic field vs temperature diagram obtained below 2 K for a:=0, 0.22 and 0.36. Solid lines 
denote the fits by eq. Q for T < 2 K. 



respectively, |27h29| . Because these two intradimer inter- 
actions are different, the partial K"*" ion substitution for 
T1+ ions produces the random local potential Vi of the 
magnon [30|. 

Single crystals of Tli_a;Ka;CuCl3 were synthesized from 
a melt comprising a mixture of TlCuCla and KCuCla in 
the ratio of 1 — a; to x. The potassium concentration 
X was determined by inductively coupled plasma opti- 
cal emission spectrochemical analysis at the Center for 
Advanced Materials Analysis, Tokyo Institute of Tech- 
nology. The specific heat was measured down to 0.36 
K in magnetic fields up to 9 T using a physical prop- 
erty measurement system (Quantum Design PPMS) by 
the relaxation method. High-frequency, high-field ESR 
measurements were performed using the terahertz elec- 
tron spin resonance apparatus (TESRA-IMR) [3l| at the 
Institute for Material Research, Tohoku University. The 
temperature of the sample was lowered to 0.55 K using 
liquid ■^He. A magnetic field was applied using a multi- 
layer pulse magnet. In all the present experiments, the 
magnetic field was applied perpendicular to the cleav- 
age (1,0,2) plane, because the magnetic anisotropy in 
this plane is so small that the U{1) symmetry is approx- 
imately conserved [33 |. 

The specific heat obtained from the temperature scan 
exhibits a small cusplike anomaly due to magnetic or- 
dering for > 4.5 T. The anomaly becomes smaller with 
decreasing magnetic field. Thus, we performed the field 
scan of the specific heat at various temperatures below 2 
K. Some examples of the measurements for a: = 0.36 are 
shown in Fig. [ijb). The specific heat exhibits a cusplike 
anomaly, to which we assign the transition field H]S!{T). 
The transition field is well-defined within an error ±0.1 
T. The phase transition points obtained below 2 K for 




FIG. 2: Low-temperature phase boundary for x =0.36 and 
critical exponent <^ as a function of Tmax obtained by fit- 
ting eq. ((2]) to transition points between 0.36 K and Tmax for 
H±{1, 0, 2). The solid line is the fit with (j> = 0.80. 



x=0, 0.22 and 0.36 are summarized in Fig.[ljc). It is 
clear that the critical behaviors of the phase boundaries 
for x — and x^O near T = are qualitatively different. 
The phase boundary for a; = is normal to the field axis 
for T but is not for xj^O. 

The phase boundaries shown in Fig.[IJc) can be de- 
scribed by the power law of eq. ([2]). For x = 0, we obtain 
(/)= 1.53 ±0.14 using the data for 0.35 <T< 2 K. The 
solid line for x = Q in Fig.[ljc) is the fit with ^=1.53. 
This critical exponent (j) coincides with 4'bec = 3/2 pre- 
dicted by the BEG theory without disorder 0, H, 0, 122] . 
On the other hand, for x = 0.22 and 0.36, we obtain 
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4> = 1.00 ± 0.06 and 0.96 ± 0.03, respectively, using all the 
data below 2 K. The solid lines for a: = 0.22 and 0.36 in 
Fig.[ljc) are the fits with these exponents. However, for 
X 7^ 0, the critical exponent (p tends to decrease with de- 
creasing fitting window. Thus, we analyze the critical be- 
havior of the phase boundary for 2: = 0.36 in detail. We 
fit eq. in the temperature range of T;„in < T < Tmax, 
setting the lowest temperature at Tmin — 0.36 K and vary- 
ing the highest temperature Tmax from 1.87 K to 0.82 K. 
The critical exponent as a function of T,„ax is shown in 
Fig. [21 The critical exponent </> decreases systematically 
with decreasing Tmax, and = 0.58 ± 0.17 for T^ax = 0.82 
K. The solid line in Fig.[2]is the fit with = 0.80 ±0.07 
obtained for Tmax = 1.1 K. With decreasing T,„ax used for 
fitting, the critical exponent 4> shows a clear tendency to 
reach a value smaller than 1/2. This critical behavior is 
consistent with that for the BG-BEC transition discussed 
by Fisher et al. [l^ . The similar behavior for the critical 
exponent is also observed for a; = 0.22. 

As shown in Fig.[IJc), the critical field He for a; = 0.22 
and 0.36 is smaller than that for a; = 0. For the fol- 
lowing reason, this should be mainly attributed to the 
fact that the average intradimer interaction decreases 
with increasing x. Within the framework of the dimer 
mean- field approximation Q , the triplet gap is expressed 
as A= [,]{,] — 2|J|)]^/^, where J is expressed by a lin- 
ear combination of interdimer interactions @. The gap 
shrinks either when the intradimer interaction J is re- 
duced or when the interdimer interaction is enhanced. 
The temperature Tyi giving the maximum magnetic sus- 
ceptibility decreases with increasing x, as shown in Fig. 1 
in Ref. [l5|. This indicates that the average intradimer 
interaction decreases with x, because Tm is given by 
1.60Tm = J/^B vly the mean- field approximation. 

To investigate the localization of magnons, we per- 
formed ESR measurements on Tli_^Ka;CuCl3 with 
x = 0.22 and 0.36. The crystals used in the specific 
heat and ESR measurements were taken from the same 
batch. We used thick samples to gain intensities. The 
measurements were performed mainly at frequencies of 
111 and 118.6 GHz for cc = 0.22 and 0.36, respectively. 
The paramagnetic resonance fields for these frequencies 
are estimated to be 3.55 and 3.80 T, using (7 = 2.23 for 
H _L (1, 0, 2). The resonance field Hq is close to the criti- 
cal field He- A single ESR line was observed in the para- 
magnetic phase for 0.55 < T < 16 K. Although the signal 
intensity decreases monotonically down to 0.55 K, the 
intensity is sufficient, even at the lowest temperature, to 
analyze the lineshape as shown in Fig. [31 The tempera- 
ture dependence of the ESR intensity is similar to that of 
the magnetization [l^ . We also measured ESR spectrum 
at 102 GHz, but we could not obtain well-defined signal 
to analyze the lineshape. 

The lineshape of the ESR spectrum is a typical sub- 
ject in the statistical mechanics of the irreversible pro- 
cess [H,!!^]. In the paramagnetic phase, the linewidth 
is produced by the local field due to the perturbation 
that do not commute with the total spin such as the 



anisotropic exchange, dipolar interaction and the Zeeman 
interaction with nonuniform g factors. From a statistical 
point of view, it is natural to assume that the local field 
has a Gaussian distribution. Thus, when the magnons 
created on the dimers localize, the lineshape should be 
Gaussian, because the distribution of the local field deter- 
mines the lineshape of the ESR spectrum. On the other 
hand, if the magnons delocalize in the crystal, the local 
field acting on spins is rapidly averaged, so that the ESR 
spectrum is narrowed and becomes Lorentzian [ssl [3^ . 
This is called exchange narrowing. 

Raw ESR spectrum is somewhat unsymmetrical with 
respect to the resonance field Hq as shown in Fig.[3ja). 
This should be ascribed to the interference of the submil- 
limeter wave inside the sample, because the wavelength of 
the submillimeter wave used and the sample thickness are 
the same order. Since the present system is insulator, this 
unsymmetrical lineshape is not intrinsic to sample. Thus, 
we symmetrized the spectrum, averaging both sides of 
Ha, i.e., we put {I{H) + I{2Ho- H)] /2 as the inten- 
sity I{H) at H . The symmetrized spectrum has the fiat 
baseline. To analyze the lineshape of these ESR spectra, 
we plot I{Ho)/I{H) against [{H - HoU(^Hi/2)]^ , where 
AHi/2 is the line width at I{Ho)/2 [3i|. When the line- 
shape is Gaussian, I{Ho)/I{H) increases exponentially, 
while for the Lorentzian lineshape, I{Hq) / I{H) is a linear 
function of [{H — Hq) / {AHi/2)]'^ with a slope of unity, as 
shown by dashed and solid lines, respectively, in Fig.[3fb) 
and (c). The I{Ho)/I{H) vs [{H - Ho)/{AHy2)? plots 
at various temperatures for a; = 0.36 and 0.22 are shown 
in Fig.[3Ub) and (c), respectively. We can see the similar 
temperature variation in these plots. I{Hq)/ I{H) at high 
temperatures is almost linear in [{H — Hq) / {AHii2)Y 
with a slope of unity, which shows the Lorentzian line- 
shape. On the other hand, at sufficiently low temper- 
atures, I{Hq)/I{H) vs [{H-Ho)/{AHii2)f plot is be- 
tween those of Lorentzian and Gaussian. Such an ESR 
lineshape is observed in one-dimensional Heisenberg an- 
tiferromagnet, in which the time correlation function of 
the local fields acting on th e sp ins does not damp rapidly 
but has the long time tail |35| . 

In the present temperature range, the BEG phase is not 
reached even at the lowest temperature. At high temper- 
atures, thermal hopping of magnons is activated. This 
causes the rapid averaging of the local fields, which leads 
to the narrowed Lorentzian lineshape. At sufficiently low 
temperatures, the thermal effect is suppressed. Thus, 
the intermediate lineshape observed for T < 0.8 K can 
be attributed to the localization of magnons. Since the 
resonance fields for 111 and 118.6 GHz are close to the 
critical field H^, it is considered that the complete local- 
ization does not occur in the present temperature range, 
so that the lineshape is intermediate between Gaussian 
and Lorentian. 

Oosawa and Tanaka [l^ reported the results of magne- 
tization measurements on Tli_2,K2;GuGl3 with x<0.36. 
The magnetization curve that they observed at T = 1.8 
K had a finite slope even for H < He. This was not 
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FIG. 3: (a) ESR spectra for a; = 0.36 observed at 0.55 K for 118.6 GHz. Lower and upper spectra are raw and symmetrized 
ones, respectively. I{Ho)/I{H) vs [{H — Ho)/{AHif2)]'^ plots for the ESR spectra measured at various temperatures (b) for 
X — 0.36 and (c) for x — 0.22. Solid and dashed lines are the same plots for Lorentzian and Gaussian functions. 



ascribed to the finite-temperature efi^ect, because pure 
TlCuCla exhibits almost zero magnetization up to the 
critical field He- This result indicates that the mag- 
netic susceptibility x — dM/dH for H < is finite in 
the ground state. Since the magnetic susceptibility corre- 
sponds to the compressibility of the lattice boson system 
K — dp/dfj,, where /i is the chemical potential, the finite 
magnetic susceptibility for T — > means that the com- 
pressibility of the ground state is finite. In the low-field 
phase below He-, long-range magnetic ordering is absent 
in spite of the finite susceptibility. These properties for 
H <Hc are consistent with the characteristics of the BG 
phase [131 ■ From these observations and the present ESR 
results, we can deduce that the ground state for H < He 
in Tli_a;Ka:CuCl3 is the BG phase of magnons. The 
gapped MI phase for H < Hb appears to be destroyed 
in the present system. 

In conclusion, from the analysis of the phase tran- 



sition data and the lineshape of ESR spectrum in 
Tli_2;K2;CuCl3 combined with the previous result for 
magnetization measurement, we demonstrated that the 
quantum phase transition at He is the BG-BEC transi- 
tion of magnons and that the critical behavior for the 
temperature dependence of the transition field near He 
is described by the small exponent characteristic of the 
BG-BEC transition 
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